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We show that the matrix elements of integrable models computed by the Algebraic Bethe Ansatz can 
be put in direct correspondence with the Form Factors of integrable relativistic field theories. This 
happens when the S'-matrix of a Bethe Ansatz model can be regarded as a suitable non-relativistic 
limit of the S'-matrix of a field theory, and when there is a well-defined mapping between the 
Hilbert spaces and operators of the two theories. This correspondence provides an efficient method 
to compute matrix elements of Bethe Ansatz integrable models, overpassing the technical difficulties 
of their direct determination. We analyze this correspondence for the simplest example in which it 
occurs, i.e. the Quantum Non-Linear Schrodinger and the Sinn-Gordon models. 
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I. INTRODUCTION 

The aim of this paper is to show that a direct correspondence exists between the matrix elements [1] 
computed by the Algebraic Bethe Ansatz in non-relativistic integrable models (in the following simply 
referred to as Bethe Ansatz models) [2] and the Form Factors considered in relativistic integrable quantum 
field theories [3]. As shown below, the relation between the two quantities can be established along the 
lines of the recent studies on non-relativistic limit of a quantum field theory [4, 5]. 

The discovery of such a correspondence may greatly help deepen our general knowledge of integrable 
models and, in particular, shed new light on the calculation of their correlation functions. The reason is 
the following: while the direct computation of Bethe Ansatz matrix elements proves to be quite a difficult 
task (often carried out successfully only for few operators), the computation of Form Factors is instead a 
simpler type of problem. In the latter case, for instance, one can take advantage of additional constraints 
coming from the relativistic invariancc of the field theory and as a result, explicit expressions of Form 
Factors can be usually found not only for a few operators but also for a larger number of them: actually, 
the classification of all operators of a quantum field theory can be obtained in terms of the different 
solutions of the Form Factor equations [6, 7]. 

At the heart of this correspondence lies the S'-matrix, both of the relativistic field theory and the Bethe 
Ansatz model: if the S'-matrix of the latter is obtained by a suitable non-relativistic limit of the S-matrix 
of the former, then the Form Factors of the quantum field theory go smoothly to the Bethe Ansatz matrix 
elements. Obviously one has to be sure that there is also a one-to-one mapping between the Hilbert spaces 
and operators of the two theories. But if one can prove that such a mapping exists, it is then easy to 
understand why the field theory Form Factors reduce to the Bethe Ansatz matrix elements: this happens 
because the analytic properties of the Form Factors and the Bethe Ansatz matrix elements are dictated by 
the S-matrices of the corresponding theories. In the following we provide evidence for this correspondence 
by analyzing the simplest models in which it occurs: the Quantum Non-Linear Schrodinger (QNLS) model 
on one side, and the Sinh-Gordon (Sh-G) model on the other. In particular, our approach provides a 
universal method to compute matrix elements of any local operator in the QNLS model. 

II. MATRIX ELEMENTS IN THE ALGEBRAIC BETHE ANSATZ 

In this section we summarize without derivation some basic properties and results of the Algebraic 
Bethe Ansatz solution of the QNLS system associated to the Lieb-Linigcr model [8]. The interested 
reader is referred to the book [2] and references therein. 
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The Hamiltonan of the model in volume L with periodic boundary conditions is given by 

L 

Hqnls = J d.T (d x ^d x ijj + c^> WV>) > (1) 
o 

where ip(x,t) and ip'(x,t) are canonical Bose fields 

[i>(x,t),^(y,t)]=6(x-y), (2) 
and c is the coupling constant. The Fock vacuum is defined as 

#))=0, (0|V f =0. (3) 
The QNLS model can be solved via the Algebraic Bethe Ansatz (BA). The monodromy matrix reads 

and its entries act in a space consisting of states 

N 

|Ai,...,Ajv> = n B ( A ;)l >> N = 0,1,.... (5) 

j'=i 

where B(A) = 5(A) exp(— iXL/2), {A} are arbitrary complex parameters, and the pseudo-vacuum |0) 
coincides with the Fock-vacuum. Similarly, dual states can be constructed using the operators C(A) = 
C(X)exp(-iXL/2). 

The i?-matrix describes the commutation relations of the monodromy matrix entries and it satisfies 
the Yang-Baxter equations. For the QNLS model it can be written in the form 

//GM) o o o \ 

V f((i,X)J 

with 

/( M ,A)= M ~ A + iC , ff ( MjA ) = ^-. (7) 

/i — A ^ — A 

The transfer matrix r(A) = trT(A) = -A(A) +D(A) generates the complete set of the conservation laws 
of the model. The eigenstates of the transfer matrix have the form (5), however the parameters {A} are 
not arbitrary but they satisfy the system of Bethe equations 

N 

zAt L 



n^QNLs(A J ,A fe ) = l, (8) 

fc=i 

where the two-particle S'-matrix is given by 

s n y \ - f(MiM A 3 - A fc - ic 
f(Xj,X k ) Xj-Xk + ic 

The S'-matrix gives the phase factor by which the state gets multiplied when the particles i and j 
are interchanged. Hence the Bethe equations say that the total phase-shift acquired when a particle of 
momentum Xj is taken to a round trip comes from the usual phase which is proportional to the momentum 
plus the scattering phase shifts picked up when the particle scatters through all the other particles. Taking 
the logarithm leads us to 



N 1 

= XjL + J log ^qnls (Xj , A fc ) = 2irlj , Ij € Z . (10) 
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Using the algebra satisfied by the monodromy matrix, the scalar products of the BA states (5) can be 
worked out explicitly, as well as the action of the operator ip on these states. However, the calculation 
of the scalar products proved to be a highly non-trivial combinatorial problem (see [2] and references 
therein). As a result, the norms of the states with parameters A that satisfy the Bethe equations (8) are 



iV 



(Ai,...,Ajv|Ai,...,Ajv) = c w [ /(A 3 -,A fc ) \ p N , (11) 




where 

is the Gaudin determinant associated to the Bethe equations (8). Knowing the action of tp on the BA 
states and the scalar products, its unnormalized matrix elements 

F$(X[, • • ■ , A'^.ilAi, . . . , X N ) = (K, . . . , X' N _ 1 \ip(0, 0)|A 1; . . . , X N ) (13) 

can be given explicitly. These matrix elements are originally defined for states which solve the Bethe 
equations. However, one can define a function Fn such that the actual matrix elements will be given by 
the value of this function taken at the particular set of momenta which satisfy the Bethe equations. Hence, 
the function Fn itself does not carry any information about the system size L: in fact, this only enters 
the Bethe equations satisfied by the physical momenta. Note that the only non-zero matrix elements of ip 
are for states where the difference of the particle numbers is one and that the functions Fn are symmetric 
separately in the momenta A and A'. We will give explicit examples for matrix elements in section IV. 
However, it is important to note that they satisfy the recursion relation [9] 

F%(X' 1 ,...,X' N -i\Xi,...,X N ) — 

(AT-1 N N-l N 

n n a i) - n fwx) n k>*> a 
3=2 3=2 j=2 j=2 

x F*_ 1 (Ai,...,A' Ar _ 2 |Ai,...,A JV _i) + ... , (14) 

where the dots stand for non-singular parts. 

Similarly, matrix elements of the current operator j(x) = ^ (x)i(j(x) can also be determined. The 
current has non-zero matrix elements only between states having the same number of particles, 

F 3 N (X[, ...,X' N \X 1 ,...,X N ) = (X' 1 ,..., A5v|V + (0, 0)^(0, 0)|Ai, ...,X N ), (15) 

and they also satisfy the recursion relations (14) with the obvious change in the number of particles of 
the dual vector (JV — 1) -> N. 

III. FORM FACTORS IN THE SINH-GORDON MODEL 

The Sh-G model is an integrable relativistic field theory in 1 + 1 dimensions defined by the Lagrangian 
density 

£ - H&)' - Ks)' <i6 » 

where <j> = <f>(x, t) is a real scalar field, mo is a mass scale and q is the speed of light. The parameter mo 
is related to the physical (renormalized) mass m of the particle by [10] 

2 2 na 

m = m — — r . (17) 
sm(7ra) 
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The intcgrability of the Sh-G model implies the absence of particle production processes and its n-particle 
scattering amplitudes are purely elastic. Moreover, they factorize into n(n — l)/2 two-body S-matrices 
which can be determined exactly via the S-matrix bootstrap [11]. The energy E(9) and the momentum 
P(6) of a particle can be written as E(9) = Mcf cosh 9, P{9) = Mc;sinh0, where 9 is the rapidity. In 
terms of the rapidities the exact two-body S-matrix is given by [12] 

b sh - G {Vi,9 2 ) = b{9 1 ,9 2 ) - . . r, (18) 

smnt'i2 + 1 sxsx[cm) 

where 9\ 2 = 9\ — 9 2 and a is the renormalized coupling constant 

hct g 2 



Sir + hci g 2 



(19) 



The key observation, made in papers [4, 5], is that the QNLS model can be regarded as a suitable 
non-relativistic limit of the Sh-G model, under which the two-particle S-matrices, the Hamiltonian and 
the Thermodynamic Bethe Ansatz equations of the Sh-G model go into the corresponding quantities of 
the QNLS model. The connection between the two theories is realized by taking a double limit, in which 
the speed of light q goes to infinity, the coupling g goes to zero, but with their product kept fixed and 
given by 

4\/c 

ci -> oo , g -> , gci= , (20) 

where c is the coupling constant of the QNLS model. Taking such a double limit of the S-matrix of the 
Sh-G model one arrives at the S'-matrix (9) of the QNLS model, once we set m = 1/2 and H = 1. Note 
that itiq — > m in the limit. 

Such a correspondence between S'-matrices gives a hint that an exact mapping exists between the two 
models. Given that in a relativistic intcgrable model the two-particle S-matrix governs its entire dynamics 
(its thermodynamic properties and the Form Factors of all its operators), if this mapping exists then it 
has several interesting and far-reaching consequences. For instance, in the past a mapping between two 
relativistic S-matrices was used to establish the correspondence between the form factors and, based on 
this, between the operators of the theories [13]. As discussed below, in the present case the situation is 
more subtle, because the mapping operates between a relativistic theory and a non-relativistic model. 
The analysis in this case requires additional care about the structure of states and operators of the two 
Hilbert spaces. 

Here we would like to shed light on the direct relation between the Bethe Ansatz matrix elements 
of the QNLS and the Form Factors of the Sh-G model. In a relativistic field theory defined in infinite 
volume the elementary Form Factors of a local operator O are the matrix elements of 0(0,0) between 
the vacuum and a set of n-particle asymptotic states: 

F^(9 1 ,9 2 , . . . , 9 n ) = (0|O(0, 0)10!, . . . ,6 n ). m . (21) 

The knowledge of all Form Factors of an operator is enough for realizing how it acts on any state 
of the theory. In fact, a generic matrix element of the operator can be expressed in terms of its Form 
Factors by using the translation operator and the crossing symmetry which is implemented by an analytic 
continuation in the rapidity variables [14]: 

(9[, . . . , 6' n \O(0, 0)^,. ..,**) = F° +k (6[ + iV, . . . , 6> n + iTr, 0i, . . . , 9 k ) . (22) 

The Form Factors satisfy a set of functional and recursive equations, which for integrable models makes 
it possible to find in many cases their explicit expressions (for a review, see [3]). For a scalar operator 
unitary and crossing symmetry dictate the following functional equations 

F n (9i, . . . , 6k, Ok+i, ■ ■ ■ , 9 n ) = S(9k — 9k+i) F n (9\, . . . , 6k+i,0k, ■ • ■ , 6 n ) , (23a) 
F n (0i + 2iri,...,0 n ) = F n (9a,...,6 n ,9i) . (23b) 

The Form Factors of integrable theories can have two kinds of simple poles and, except for these singular- 
ities, they are analytic in the strip < Imfry < 2ir (here 0y- = 9i — 9j). The first kind of poles corresponds 
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to kincmatical singularities at 9ij = in and their residues give rise to a set of recursive equations between 
the n-particle and the (n + 2)-particle Form Factors 

F n+2 (9' + i7r,e,e 1 ,...,e n ) fi-JJ 5(0-0,) ] F n (e 1 ,...,e n ) + ... , (24) 

where the dots stand for regular parts. The second kind of poles is related to the bound states of the 
theory, but since there are no bound states in the Sh-G model, there are no such poles in the Form Factors 
of this theory and we will not consider them any further. 

Note the striking similarity between the recursive equations (24) and (14): we will show below that, 
indeed, they exactly correspond one to the other in the limit (20). These recursive equations, together 
with the requirement of the correct asymptotic behavior and of the desired analyticity properties, are 
the key tools in finding explicit solutions for the Form Factors. In the Sh-G model a concise expression 
is provided by the Form Factor of the exponential operator [15] 

F n (k) = <0|e**|*, ...,9 n ) = ^« ? det M n (k) f[ . (25) 

sm(7ro;) V F m i n (iir) J e"* + e a ' 

j<i 

Here k is an arbitrary real number, F m i n (9) is the minimal solution of the Form Factor bootstrap equations 
and M n is a (n — 1) x (n — 1) matrix 



(n) sin((j-Z + fc)7ra) 
[Mn{k)] ^ = ^a) ' (26) 



where Oj are the elementary symmetric polynomials of the variables e° 



of 



J2 eV..e*». (27) 



We will see that the Form Factors of e kg ^ act as generating functions for the Form Factors of all the 
powers of the field (j>. 



IV. QNLS MATRIX ELEMENTS FROM SH-G FORM FACTORS 

Before going on with the analysis, it is worth emphasizing the strong similarities of the QNLS and 
Sh-G models, in particular, the similarity of the key equations for the Bethe Ansatz matrix elements 
and the Form Factors. Both theories are integrable and they contain a single type of massive particle 
without additional bound states. The pseudo-vacuum, on which the BA states are built, coincides with 
the Fock-vacuum, i.e. the zero-particle state. The Sh-G Form Factors are built on the physical vacuum, 
but in the non-relativistic limit the zero-point fluctuations disappear and we obtain the Fock-vacuum of 
the QNLS model. The two-particle S'-matrix, which in the relativistic context governs every aspect of 
integrable theories, also corresponds between the two models. From this it is clear that the Bethe equations 
(10), (29) are also mapped to each other. The striking similarity between the recursive equations (14) and 
(24) is very important, because these are the key equations that allow for the determination of the Form 
Factors and the matrix elements. Moreover, the building block of the relativistic Form Factors, F m i n (9) 
has the following behavior: in the limit (20) 

iWfy) -> i A " 7 A ' ■ = /(A,-, AO , F mia (6 jt for 9 n g R , (28) 

Aj — A; + IC 

thus it goes into an important building block of the Bethe Ansatz quantities. With all these correspon- 
dences it should be quite clear that the Sh-G Form Factors will be mapped to the QNLS matrix elements. 
In what follows we give the details of this mapping and provide explicit examples. 

In order to obtain QNLS matrix elements from the Sh-G Form Factors using the limit (20), we have 
to understand both the relation between the operators and the states of the two models. We also have 
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to take into account that the Bcthc Ansatz matrix elements are defined in finite volume whereas the 
relativistic ones discussed so far are defined in infinite volume. The latter difference can be cured using 
the results of [16]. Even in relativistic quantum field theories the states in a finite volume L can be 
described as multi-particle scattering states, where the rapidities of the particles are solutions of the 
Bethe quantization conditions: 

n 1 

Q 3 = mciL sinh *,■ + ^ - log S(9j - 9 k ) = 2nlj , j = 1, . . . , n . (29) 

However, this procedure is to be understood as an "asymptotic Bethe Ansatz" in the sense that it gives 
the correct results to all orders in 1 / L but it misses residual finite size effects which decay exponentially 
with the volume. These corrections can be associated to processes involving virtual particles with some of 
them "travelling around the world" . The main idea of [16] is that this picture of "asymptotic BA" applies 
to the form factors as well: the matrix elements in a finite volume are given by the infinite volume Form 
Factors at the particular set of rapidities which solve the corresponding Bcthc equations. In addition, one 
has to introduce normalization factors given by the corresponding Gaudin determinants 

p n (e 1 ,...,e n ) = det^- (30) 

which can be interpreted as the density of states (in rapidity-space) in the corresponding sector, or 
alternatively as the norms of the BA states like in (11). Thus the finite volume Form Factors are given 
by the infinite volume ones taken at the rapidities satisfying the Bethe equations (29), divided by the 
density of states: 



<*;,..., 310(0, OH*, ... , e n ) L = f°^ + ^---A + ^...^ ) L 



where 0(e~ flL ) stands for the above-mentioned exponentially small corrections in L. In the non-relativistic 
limit both the Bethe equations and the norms of states go over to the QNLS quantities (10), (12). In 
particular, the relation between the norms is given by 

Pn — > (mci) N p N . (32) 

We note that (31) is valid as long as there are no coinciding rapidities in the "bra" and "ket" vectors. In 
the latter case disconnected terms arise; a proper treatment of such contributions was given in [17]. In 
this paper we do not elaborate on diconnected pieces, i.e. we assume that the two sets of rapidities are 
completely distinct. 

Let us turn now to the relation between the operators in the two theories which was given in [4]: 



<f>(x,t)~ \ — (^OM)e-^ *+^ f {x 1 t)e +l ^ t \ . (33) 
V 2m \ / 



The exponential terms have to be separated, because the relativistic Hamiltonian contains also the rest 
energy which is absent in the non-relativistic case. The sign ~ means that in any functional expression of 
4> the surviving exponential terms should be dropped, because in the non-relativistic limit c/ — > oo they 
are rapidly oscillating and give zero when integrated over any small but finite time interval. 

Similarly, we must compensate for the rest energy in the time evolution of the states, so (without the 
proper normalizations) 

|*i,. ..,*„) <— > e- inmc ' t \X 1 ,...,X n ) , (*i,...,*„| <— > e +mmc ' t (A 1 ,...,A„| . (34) 

The relation between the rapidities and the momenta in the non-relativistic limit is Xi = 6i/mci. Note 
that the encounter of the oscillating terms in the states (38) and in the operator (33) (as well as in all 
its powers) will ensure that in the non-relativistic limit a given operator (e.g. ip) have non-zero matrix 
elements only between states in which the number of particles differ by a fixed amount (e.g. one). 

Let us now deal with the question of the normalization of the states. First of all, one should note that 
in contrast with the BA states the relativistic asymptotic states are not symmetric in the rapidities, they 
rather obey 

1*1, ■■■,9k: 9k+l, • ■ • , 9 n ) = S(8k — *fe+l) |*l, ... , *fe+l, 9k, ■ • ■ , * n ) • (35) 
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Hence, in order to establish the correspondence with the BA states, these states should be symmetrized 
in the rapidities, which can be done by multiplying with the appropriate phase factors: 



'1 1 • • • i u n, 



symm 



lJs(e j -e k )\e 1 ,...,e n ) . (36) 



j>k 



The normalization (11) of the BA states should also be taken into account. From (32) it is clear that for 
the proper normalization we have to include a factor of y'rac; for every particle. 
Collecting everything we finally arrive at the relation 

J] s js{9 j -9 k )\9 l ,...,9 n ) { ™ Cl) n/2 e"hAi,-..,A w ) . (37) 

J",fe=l 

Taking the S-matriccs to the right hand side and using eqn. (9) for the S'-matrix recovered in the double 
limit of the Sh-G model, this can be written as 

(mci) n / 2 e~ in mc ^ t 

\9i,...,6 n ) ~ — |Ai, . . . , A„) , (38) 

c" /2 I! /(A 3 -,A fc ) 

j<k 

The relation between the dual vectors is given by the complex conjugate expression, in particular, the 
sign of the time-dependent phase (34) is opposite. 

Using the ^-matrix (9) of the QNLS model and the correspondence of the states (38) it is easy to 
prove, as shown below, that the relativistic recursive equations (24) transform exactly into the recursive 
equations (14). Note also that due to (28) the /(Ai,Aj) factors in (38) will exactly cancel the limiting 
forms of F min (9ij). 

Now we are in a position to obtain a generic QNLS matrix element of the form 

(A' 1 ,...,A^ +p _ ? |V tj Wi,...,A Ar ) (39) 

by performing the following steps: 

1. We determine the infinite volume Sh-G Form Factor 

(0| : r +q :\9' 1 + i7r,..., 9' N+p _ q +in,9 1 ,..., 9 N ) 

by picking the 0(k p+g ) term in the expansion of (25) and (for p + q > 2) by taking into account the 
normal ordering issues [18]. The pre-factor in (33) is 1 with the choice m = 1/2, ft = 1, but simple 
combinatorial factors may arise. As we said, the number of crossed rapidities will automatically 
select the correct combination ip^ p tjj q out of <f)' p+q . 

2. We calculate the finite volume Form Factors using (31). However, in the double limit the norms 
of the finite volume states will go to those of the BA states (11) up to a factor of ^Jmci for every 
particle, which needs to included. 

3. We take the double scaling limit (20) using (19) as well as the relation 6*,; — > A;/toq and the limit 
of the minimal Form Factor F m [ n (28). 

4. We take into account the different normalizations according to (38). 

5. To get a proper matching with the Bethe Ansatz matrix elements we include a factor of —i for each 
ip and a +i for each ip'. 

These steps can be summarized in the formula 
(Ai,...,A^ +p _ 9 |^t^| Al) ... )Ajv ) = 

/ I \ - 1 /O \ 4* N N+p-q 

'"(":") I? " .^HftA) n /«.*>* 

y J v J j<k j<k 

x \Tm{(mc l )-^^(Q\:r +q -A0[+™,... 1 6' N+p _ q + tn,9 1 ,...,9 N )} , (40) 
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where lim denotes the double scaling limit (20). 

As a first check of this procedure let us explicitly show the correspondence of the recursive equations 
(24) and (14). Starting from the relativistic case we have 



F 2N -i(0'i +in,...,9' N _! +in,0i, 



?7VJ = 



N-l 



I S{9' 3 + ITT - 9 1 )F 2N . 1 (9[ + ITT, 9 U 9' 2 + tTT, . . . , 0^-1 + ITT, 2 , . . . , 9 N ) 



3=2 

AT-1 AT 1 N-l 



1- II S{6' l -9' j -m)J[S{9 1 -9 j ) f II S(^^-Wjv-3(^+*'r.-.- ) ^-i+»'r^a I ... J M 

j=2 j=2 



JV-1 A? 

*i - 6»i 

In the limit using (40) this relation becomes 



H S(9[-9' j )-l[S{e 1 -6 j )\F 2N _ 3 {6' 2 + iTr,...,9' N _ 1 + iTr,9 2 ,...,9 N ) . (41) 

3=2 j=2 



(^) 2 V^"^ - ~V — ^ A i""' A ^-i' Al '---' Ajv ) 

n /(A 3 -,A fc ) n /(^A' fe ) 



a; - Xi 



AT-1 N 

J -5qnls(A' 1 — Aj-) - S*qnls(Ai - Aj) 

3=2 3=2 



x ? i 1 ?) 2Nf3 v£— — ^ — ^- i(a - - aw - i; A2 ' ■ • ■ - xn) • (42) 

n /(Aj,A fc ) n /(a;-. Ay 

2<j<k 2<j<k 

Simplifying with the common pre-factors and using relation (9) between the S'-matrix and the function 
/(A,ju) we arrive at equation (14). 

Let us take now some explicit examples of the limit for the matrix elements of the field operator -0 and 
the current operator ip'ip. For the matrix elements of ip a nice determinant representation was found in 
[19]. They can be expressed also as [9] 

pVvy y u \ ^ _ fiv(Ai,...,A / jV _ 1 |Ai,...,AAr) 

-f , A r(A 1 ,...,A Ar _ 1 |Ai,...,A A r) - — , {46) 

n n(A J -A' fc ) 

fe=l j=i 

where Pjv are polynomials in {A}. The first few examples are given by 

Pi = -iy/c , (44a) 

P 2 = -2i^Zc 2 , (44b) 

P 3 = -4*v^c 4 [-c 2 + ((Ai - Ai)(A a - X' 2 ) + (Ai - A' 2 )(A 3 - A' x ) + (A 2 - Ai)(A 3 - A' 2 ))] . (44c) 
The form factors for ip^ip can be extracted from the matrix elements of the non-local operator 

X 

Q x {x) = J dyj(y) (45) 
o 

by differentiating with respect to x. Matrix elements of Q\{x) are listed for example in [20] and a 
determinant formula can be found in [21], which yield 

F/(Ai|A 1 ) = c, (46a) 

HVr^XyM = -2c 3 (Al + A2 ~ A,1 ~^ )2 , (46b) 

11 (A; - A fc ) 

3,k=l 
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and so on. 

On the Sh-G side the Form Factors of the various powers of <f> can be obtained by series expanding 
formula (25) in the auxiliary real variable k [18]. For example, the Form Factors of <f> are given by 

(CM*, ... , n ) = ^ (jr^) * (sin^a))-- 1 det M n (0) f[ , (47) 



or explicitly 



3<l 



mdl) = iv c \ /■ t \ ' (48) 

3 



h,e 2 ,e 3 ) = _ * ™vSR^ n ~^T~r~9T^ ( 49 ) 

Now we can apply the rule (40) with p = 0, q = 1 and iV = 1, 2, . . . It is useful to note that 



2c na y/c 

Cl a ~T ' o • 50 ) 

TTh g 2n 

From (28) we see that F m - m (iTv) — > 1 and that the surviving F m - m factors exactly cancel the /-functions 
appearing in (40). Taking care of the pre-factors the double limit yields a final result which coincides 
exactly with the expressions (44). 

Similarly, the first Form Factors of <fi 2 obtained from (25) are 

W 2 |gi,g 2 ) = F V ; /'f r iWfli-ft) - (51a) 
^min(«7r) 5 Z sm(7ra) 

(01^1^,^,03^4) = -i^^n ^.^-^ X (51b) 

^min(«7r)^ g l fJ- e w J + e w < 

x ( e fli+«»+«s+»*( e * + e ^ + e ^ + e ^)2 _ ga^+fe+fla+sO^-fli + e -e 2 + e -e 3 + £ -e^ (51c) 

Applying the double limit formula (40) with p = q = 1 to (51) one retrieves the matrix elements (46). 
We have also checked the matrix elements with higher number of particles finding a perfect agreement 
with the corresponding Bethc Ansatz matrix elements. 



V. CONCLUSIONS 



In this paper we have analyzed the close correspondence between matrix elements of Bethe Ansatz 
models and Form Factors of relativistic integrable field theories. If the former models can be regarded 
as non-relativistic limits of the latter theories, then the Bethe Ansatz matrix elements can be efficiently 
obtained as the non-relativistic expressions of the corresponding Form Factors, whose explicit computation 
is much simpler. 

We have discussed in detail this correspondence between the Quantum Non-Linear Schrodingcr model 
and the Sinh-Gordon model, where we gave a universal method to compute all the matrix elements of 
every local operator, but there are strong arguments in favor of its validity for other pairs of models 
as well. Indeed, both in Bethe Ansatz models and quantum field theories the main properties of matrix 
elements are dictated by their S'-matrices: therefore, if there is a mapping between the Hilbert spaces and 
operators of the two theories and the expressions of the two S-matrices coincide in the non-relativistic 
limit of the quantum field theory, these two facts induce a mapping between the matrix elements of the 
two theories. The use of Form Factors may help in solving most of the technical obstacles that have 
prevented so far the computation of matrix elements in non-relativistic Bethe Ansatz integrable models, 
thus opening new perspectives on the computation of their correlation functions. Along this direction it 
would be interesting to investigate multi-component systems, where the BA results are limited due to 
the nested nature of Bethe Ansatz. 
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